Symmetry adapted finite-cluster solver for quantum Heisenberg model in 
two-dimensions: a real-space renormalization approach 



C3 



I 

T3 



V.E. Sinitsyn, I.G. Bostrem, A.S. Ovchinnikov 

Department of Physics, Ural State University, 620083, Ekaterinburg, Russia 
(Dated: February 6, 2008) 

We present a quantum cluster solver for spin-S Heisenberg model on a two-dimensional lattice. 
. The formalism is based on the real-space renormalization procedure and uses the lattice point group- 

• theoretical analysis and nonabelian SU(2) spin symmetry technique. The exact diagonalization 

' procedure is used twice at each renormalization group step. The method is applied to the spin- 

half antiferromagnet on a square lattice and a calculation of local observables is demonstrated. A 
symmetry based truncation procedure is suggested and verified numerically. 

I. INTRODUCTION 

3' 

Low-dimensional magnetic systems are currently a subject of intensive experimental and theoretical work. Cluster 
methods, which approximate the physics of the infinite system by solving the problem for a corresponding finite 
cluster, are the most frequently used theoretical approaches, as they account short-range correlations on the scale of 
\ the cluster size. Numerical standard methods in the field, such as quantum Monte Carlo (QMC), exact diagonalization 
(ED)i, and density matrix renormalization group (DMRG^i 3 -, are able to give essentially exact results on limited size 
c/3 I systems and form a versatile methodological triad in simulations of model Hamiltonians. 

Even though these techniques have had spectacular successes in calculating ground state energies and many other 
properties of one-dimensional (ID) and two-dimensional (2D) quantum spin systems 4 "" 5 -"^ there is a problem with an 
utilizing of symmetries and good quantum numbers of the Hamiltonian, which may be exploited to thin out Hilbert 
space by decomposing it into a sum of sectors. Common symmetries and conservation laws encountered in spin 
systems are: (i) Ising or XY symmetry (magnetization conservation 5 t z ot =const); (ii) point group symmetry (parity, 
angular momentum conserved); (iii) full SU{2) symmetry (5 2 ot conserved). Among these symmetries only the first 
is usually exploited in numerical calculations. The full 5(7(2) spin symmetry is rather hard to implement, since 
it requires efforts similar to the diagonalization of the actual Hamiltonian to construct the eigenstates of 5 2 ot . An 
' implementation of nonabelian SU(2) spin symmetry based on Clebsch-Gordan transformations and elimination of 
. quantum numbers via the Wigner-Eckart theorem was performed for the interaction round a face (IRE) models in 
the framework of the IRF-DMRG method^. This technique has been successfully applied to the spin- 1/2 Heisenberg 
chain and, later, to the spin 1 and 2 Heisenberg chains 9 . The performant DMRG method conserving a total spin 
quantum number has been suggested by McCulloch and Gulasc i 10 ^ 11 . An application of SU(2) symmetries for the 
matrix product method (MPM) closely related to the DMR G 12 ' 13 gives a rotationally invariant formulation valid for 
■ spin chains and ladder a 13 ' 14 . 

As for the lattice point symmetry, despite its importance in characteristing energy states of a spin system, there 
I""-; , appears to have been little previous work on the subject. Even though an implementing this symmetry does not lead 
to a drastic reduction of a dimension of the Hilbert space sector to be diagonalized we can resolve properties as a 
function of additional quantum numbers (irreducible representations of the point group). This circumstance might 
be crucial for efficient truncation of Hilbert space in algorithms based on real-space renormalization group (RSRG) 
procedure. This generates a motivation for the present paper, namely, we present a finite cluster solver based on 
RSRG scheme which allows to exploit both the continuous nonabelian SU (2) symmetry and discrete symmetry of the 
lattice point group in application to isotropic two-dimensional spin-S systems. As an example illustrating features of 
our method we consider the spin-1/2 Heisenberg antiferromagnet (AFH) on a square lattice. This choice is motivated 
by two reasons. First, the physical properties of the 5 = 1/2 AFH model on the infinite square lattice at T = 
have been much studied and calculated by various methods by many physicists 15 '- 1 ^. The focus has generally been 
on the ground state energy and staggered magnetization although some other quantities have also been computed 
(see Refjii, for example). Second, the underlying idea of our approach was first developed by Lin and Campbell in 
the study of this model syste m 18 i 24 . Before moving on to the details, we discuss important aspects of finite-lattice 
simulations using the ED method regarding the cluster geometry. 

The method of exact diagonalization has been used on the best bipartite finite square lattices with up to N = 38 
vertice o 25 i 27 . On each of a set of finite lattices, the Hamiltonian of the quantum spin model is diagonalized exactly to 
find the ground-state energy and the ground state eigenvector. The ground state properties can then be calculated 
exactly. The exact ground-state data for each physical property of the model on all finite lattices are extrapolated 
against an appropriate inverse power of N to obtain an estimate of the property on the infinite lattice at zero 
temperature 2 ^. Haan et ali 3 ^ showed that certain parallelogram clusters could produce good results in finite-size 
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exact diagonalization calculations. Later, Betts et al. developed a grading scheme of parallelogram tiles of the square 
lattice that could generate the best finite clusters 2 ^. From a symmetry point of view, this approach has an apparent 
flaw: the point symmetry of a parallelogram cluster does not match that of the infinte square lattice. 

In this respect, a renormalization-group (RG) approach suggested by Lin and Campbell combining exact diagonal- 
ization results with a RG type analysis seems us to be more promising. The basics idea of their calculations for 2D 
AFH on an n x n cluster (n is an odd integer) is to divide this cluster into two parts: an inner (n — 2) x (n— 2) 
cluster and the perimeter ("outer ring"). The AFH model is firstly solved for the inner cluster and its ground state 
is mapped onto a single effective spin (all the excited states are thrown away). Thus, the problem is reduced to an 
effective ID AFH model in which spins on the outer ring experience antiferromagnetic nearest-neighbor interactions 
and interact individually with the effective central spin. The latter plays the role of a staggered external magnetic 
field. The procedure is repeated for increasing values of n and demonstrates convincingly that the staggered magnetic 
long-range order exists at zero temperature. We note especially that the lattice point symmetry holds for all the 
clusters and their ground state always has spin S — 1/2 (not a singlet) according to the Lieb-Mattis theorem^!. 

We offer to change real-space RG strategy of the approach making it closer to DMRG methodology Our treatment 
begins by dividing a cluster into a central spin and its environment. In the course of real-space RG iterations the 
environment increases (technical details are discussed in the text) and we determine how coupling between the central 
spin and the environment varies. Note especially that we address to the exact diagonalization procedure twice at each 
RG step. The first use gives access to states of the environment and the second one does to those of the whole cluster 
that provides its spectrum and observables of interest. Within our RG framework, local results such as the energy 
per bond e and the staggered magnetic moment m are measured on the central site. 

We have carried out the renormalization procedure through systems of size y/l7 x y/17, and, in contrast to approach 
in RefJ^, we keep not only the ground state. For small clusters (y/E x y/E, 3x3, y/l3 x V13) we use all of the 
excited states of the environment found by exact diagonalization, whereas for the cluster of size y/Vf x y/Vf we apply 
a symmetry based truncation procedure, retaining only the states with largest weight in the environment density 
matrix. For this cluster we have compared the exact diagonalization result for e and m with those obtained via our 
renormalization group, and we regard the resulting better than 10 -2 % agreement as support for the reliability of our 
calculations. 

We note that several other methods to improve the RSRG calculations have been previously formulated to study 
low-energy properties of spin lattice models. Among the most important and successful ones, one may cite the Real 
Space Renormalization Group with Effective Interactions (RSRG-EI}i2, and the Dressed Cluster Method (PCM) 20 ' 21 . 
The first method is an improvement of the RSRG method originally proposed by Wilson. By considering the blocks of 
lattice it extracts effective interactions between the blocks through the exact diagonalization of dimers of blocks. The 
knowledge of the exact spectrum of the dimers enable one to define interblock effective interactions via an effective 
Hamiltonian. This procedure is iteratevely repeated to blocks of blocks providing at very low cost reasonable estimate 
of the energy per bond for ID and 2D spin lattices. The second method (DCM) uses a single reference wave function 
as do the Coupled Cluster Method (CCM) 22 i 23 . This wave function is used as a bath in which a finite cluster is 
embedded and treated exactly. The effect of excitations occuring on the bonds around the cluster is taken into 
account through a dressing of the cluster configuration interaction (CI) matrix. This approach gives results for the 
cohesive energy of the same accuracy as the best QMC ones. The DCM can be seen as a convenient approximation 
of the CCM. However, the problem is formulated as a diagonalization of a dressed CI matrix instead of the resolution 
of a nonlinear system of equations. 

The paper is organized as follows: The general formalism for two-dimensional spin-S" systems is introduced in the 
following section. In Sec. Ill we apply the method to the AFH model on a square lattice. Finally, our conclusions 
and an outlook are presented in Sec. IV. 

II. TWO-DIMENSIONAL ISOTROPIC HEISENBERG SPIN-s SYSTEM. 
A. Cluster states and observables 

In the first step one must identify the cluster. As detailed above, care should be taken to ensure that the cluster 
has the same point-group symmetry as the lattice. Since calculation of antiferromagnetism requires bipartite clusters, 
we select a cluster with a bipartite environment of the central site (the case of this violation will be illustrated in the 
example of the cluster y/l3 x VT3). 

The cluster Hamiltonian 
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(i) 



is composed of the term V = J So ^2 Sq+s describing interactions of the central spin So with the nearest neighbors at 

distances S and rest terms denoted as the Hamiltonian of the "environment" H u . Since, by construction, the cluster 
retains a lattice point symmetry, its states \iSMTfi) with the energies Eisv are labeled by the cluster total spin S 
with the third component M and by the irreducible representation Tfj, of the cluster point group. Different states 
with the same values SM and Tfi are distinguished by the index i. In addition we need to consider the operator 

Oqi 1 = S (^Q+i) as a double irreducible tensor which transforms according to identity representation A\. The 
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same arguments enable us to use the irreducible form of the central spin operator (So)* = (So)qf 1 ■ The part V may 
be written as the inner product 

v = jJ2(-i) 
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where 



M Ax Ax 
1 1 1 



= 1 is the Clebsch-Gordan coefficient of the cluster point group . 



Let us suppose that we have found the eigenvalues -E^s^r,, and the eigenstates of the environment Hamiltonian H u 
in the form \i u S u M u T u fj, u ), The basis functions of the full cluster are obtained by the addition rule of spin angular 
momentum 



\iuS u F u ] s\ SMT u /j u ) — 



flu ,cr 



S u s S 
M u a M 



\i u S u M u T u fj, u ) |ser) 



(2) 



where [. . .] is a Clebsch-Gordan coefficient, hereinafter we use that of given in Ref^, and \sa) is the wave function of 
the central spin. Since the state \sa) is invariant under all transformations of the point symmetry group, the cluster 
basis functions transform like that of the environment according to the same irreducible representations. 

The calculation of matrix elements for the Hamiltonian (TT|) with the help of the Wigner-Ekart's theorem yields (see 
Appendix A) 

i u S u T u ;s;SMT u ii u H i' u S' u T' u ; s; S' M'T' u /j,' u 



E, 



i„s„r 



A^'Js u ,s>Jt u t>J^>Js,S'5m,m' + J(-l) s:+5+1/2 



S u s S 
s SL 1 



M,M' 



x (s \\s\\ s) (i u S u T u \\O lM || i' u S' u T' u ) Sr^r'J^, 



(3) 



where {...} is a 6j-symbol. The first reduced matrix element is (s ||s|| s) = y^s(s + l)(2s + 1) and the latter may be 
obtained if the environment eigenstates are known (see Sec. III). The energy per bond is then calculated as 



£ isr„ = 



z V 6c 



= - E 



E 



(env) 



where z is the number of nearest-neighbors of the central spin. The eigenfunctions 

\iSMTfj,) = PiOs u r u \iuSuT; s; SMTfj) , (Tfi = T uf i u ) 



(4) 



(5) 



and the energy levels -Eisr are determined by direct diagonalization of the cluster Hamiltonian H [Eq.@]. The values 
Eisv should be regarded as an approximation of the energy spectrum in the thermodynamical limit, whereas the 
energy Eisr divided per bond number is much less appropriate for this. 
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It is important to note that from Eq.Q it follows that to build the cluster target state \iSMT/j,) we need only to 
know the states \i u S u M u T u fi u ) of the environment with the quantum numbers \S — s\ < S u < S + s and T u fi u = Tfi. 

The most important quantity typically measured in numerical simulations is the ground-state staggered magneti- 
zation M c . The quantum mechanical observable for z-projection of the central spin is given as follows 



(iSMTfi \S%\iSMTn) = (-l) 1+s+s Mi 



/ 25 + 1 
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where the identity (|32|) is used. The staggered magnetization M c is determined as 



M c 2 = lim 3 



S z (R)S z (0) 



where factor 3 arises from rotational symmetry in spin space. At long distances (^S z (R)S Z (0)j ~ (5 Z (0)) that 

yields our estimate of the full root-mean-square staggered magnetization per spin M c — *J 3 (Sq) 2 . 

According to Eq. (|25p . spin-correlation function in the states of Ax -symmetry, the ground state symmetry as shown 
below, is determined as 



(iSMA 1 \S^S Z \ iSMA^ = - (iSMAi 



SoSj 



iSMAi 



= 3^ E E €sU^tstM (- 11*11 {iuSu \\S' A (rjJH i' u S' a ) (-1)*+^ { ^ £ f X , (7) 

where z/ is the lattice coordination number. In this calculation it is convenient to introduce the double irreducible 
tensor S^i (rj) — (Sj)^ summing spins at distance rj, which transforms according to identity representation A%. 
One can see that O 1 ^ 1 = S\f (6) . 



B. Increasing cluster size 



As mentioned above, the lattice point-group symmetry should be conserved with increasing cluster size. The 
requirement is put into a practical computational scheme by the following algorithm: (i) At step N we have the 
eigenvalues -E^s^r,, ano - eigenvectors \i u S u Tn u T u ^, u ) of the environment. Make a regular symmetry conserving 
expansion in the cluster size by adding sites from the next coordination shell, (ii) Using a scheme of coupling of angular 
momenta we build the set \iiSimi) of states with total spin Si and third component mi for the part that is being 
attached to the environment. The index ii labels other possible quantum numbers, (hi) In general case, these functions 
form a basis of reducible representation of the cluster point group. Based on the projection operator technique, one 
build basic functions \iiSimiTifii) transforming according to irreducible representations IV/i/. (iv) Using a scheme of 
coupling of angular momenta build a new set \iiSimiinSnmn; S u m u T u fj, u ) (jv +1 ) of states associated to the extended 
environment, where the notation \inSiimnT nfJ-n) = \iuS u m u T u fi u ) is introduced. An interaction between the 
7V-th step environment and the part added to it can be conveniently written through the irreducible tensors U ltul 
and W u<>1 built from spin operators of the "old" and "new" added parts, respectively, 

^EEEM) 9 

t„t q 

The indices t u to label different tensors of the same symmetry. The matrix elements of the extended (N + l)-th step 
environment is 

(iiSimiiiiSiimii; S u m u T u fi u ^^i'jS'im'ji'jjS'um'jj] S' u m' u T' u n' u ) = ^s„sv <5 mum ^<5r u r; 4 ^„^ 
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(8) 
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The derivation of Eq. |8j and the reduced matrix elements of the operators involved in Eq. ((SJ are given in Appendix 
A and Appendix B, respectively. 

At final step, we diagonalize ([SJ and find the eigenvalues E^ N J~^ and eigenvectors 

The iteration is closed by recalculating reduced matrix elements of the irreducible tensors W lfin in the basis of the 
extended environment (see Appendix B). Note that following the scheme we will in some cases form intermediate 
clusters, unsuitable for calculations of local results, with a non-bipartite environment. 



Si Sh S u 
mi mu m u 



Ti Tn T u 
Hi Hii Hu 



\iiSimiTifii) \inSnmiiTunn) 



III. AN EXAMPLE: SPIN-1/2 ANTIFERROMAGNET ON A SQUARE LATTICE. 



The spin-half antiferromagnet on a square lattice represents an optimal playground to study the strength and 
limitations of the method. To implement the algorithm, we need first to build wave functions of the environment 
which are predetermined by the lattice point symmetry. 

To perform calculations we start with the cluster of minimal size x a/5- The sequence of clusters involved in the 
calculations are shown in Fig. f . Within the smallest cluster, the central spin interacts with the nearest environment 
consisting of the spins S ai , , S 7l , S ril . The spin wave functions of the environment with the total spin number S u 
and the third component M u may be written as follows 



E E 

m ai ,m l 3 1 ,m T1 ,m ni m.^j ,m~, ini 



1/2 1/2 S aiPl 
m Ql mp 1 m aij 3 1 



1/2 1/2 S Jim 



ni-. 



x \l/2m ai ) \l/2m f3l ) |l/2m 7l ) 11/2™^) . 

In such a description, all allowed configurations are comprised by a set |00; 00), 1 1 1 ; 00) , |01; 1M), 1 10; 1M), |11; 1M), 
1 11; 2M), where we have dropped the spin 1/2 arguments for notation convenience. It is easy to see that the functions 
\S ai j3 1 S lim ; S U M U ) form (in common case) a basis of reducible representation of the group D4 (for details see Appendix 
C) 

g \S ai f3 1 S~ /1 r ll ; S U M U ) — ^ > s'^ i0i Si tin:L ,s ai01 Sy ini (9) \S ai/ 3 1 S liril ; S U M U ) . 

The matrices D^ K ,{g) (the upper index denotes the spin S u ) with the multiindex n — {S ai /3 1 S lirjl } are readily 
determined and read 



1(0) 



(0) 



,(0) 



1 
1 



,(o) 



(0) 



,(0) 



1/2 -V3/2 \ 

-V3/2 -1/2 y 



The functions |00; 00), 1 1 1 ; 00) form a basis of this two-dimensional representation. Still another representation of 
D4 can be generated by means of the functions |01; 1M), 1 10; 1M) and 1 1 1 ; 1M) 




-1/2 -1/2 1/V2 \ 
0JS'(Cl)= I -1/2 -1/2 , 

-1/V2 1/V2 / 

-1/2 -1/2 -1/V2\ 

^i!i'(^4)= ( -1/2 -1/2 1/V2 , 
1/V2 -1/V2 / 



1/2 1/2 1/V2 \ 
DZ,(C%)= I 1/2 1/2 -1/V2 

1/V2 -1/V2 / 



1/2 1/2 -1/V2' 
1/2 1/2 l/\/2 

■1/V2 1/V2 
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-1 \ / -1 

<UO= ( -10 , D^,(a' v )=[ 1 
1/ \ -1 

(2) 

In a similar way we find the matrices D^/ K (g) in the basis |11; 2M) 

D^M = 1 (Vs 6 D 4 ). 

The representations are the direct sums of the irreducible representations = D^ 0Al ^ © D^ 0B2 \ = 

£)( 1B i) q D^ 1E \ — D( 2Ai ^ (see Appendix E). The basis functions of these irreducible representations are given 
by a similarity transformation 

\S u M u ;Tfj) = ^ ^s^f^s^-Xfj. ISai^S^m; S U M U ) , (9) 

and the matrix u ' q , r —Tq _ o ^ o /-s«r«u»^5« S' found with the aid of the projection-operator technique 
reads (for details see Appendix D) 

|00;Ail) |00;B 2 1) |lM;Bil) \1M;E1) \1M;E2) \2M;Axl) 
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reduced matrix elements of the 
double irreducible tensor 1Al = S ai + Sp 1 + S 11 + S, n can be computed straightforwardly using the Wigner-Eckart 
theorem and the similarity transformation §§§ 
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(5r||0^||5T'>= E ^Ss^^r^s^s^sTV^^^^HO 1 !!!^^;^), (10) 

Si2,5 , 34 s[ 2 ,s' 3A 



where the indeces a±, 0i, 71, 771 are correspondingly denoted by the numbers 1-4. 

To calculate the reduced matrix element that comes into the right-hand side of Eq. (flT))) one has to rewrite O^, 
through the spin operators and employ their expressions for the reduced matrix elements of the spin operators 



(S12S34; S \\Si\\ S' 12 S' 34 ; S') — ^_2) 1+Sl2 + 534 + 5l2+,s [S12, S' 12 ,S, S'] 



1/2 
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1/2 1/2 1/2 J \ S 12 5 34 S[ 2 
(S 12 S u ;S\\S 2 \\S[ 2 S 3i ;S') = { ^+^+s 3i+ s' [<? 12j <^ 2 , S'] 1/2 

(l/2||5||l/2)^ 34 ^ 4 , (12) 

(S12S34', S \\S 3 \\ S' 12 S 34: ; S') = (_iy+ 5 i 2 +- s 34+ s s' 3i , s, s'] / 

(l/2||5||l/2)^ iaiSi3l (13) 



S' l2 1 S 12 \ J S' 1 S 
1/2 1/2 1/2 / 1 S 12 S 3i S' 12 



x 
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TABLE I: Energies Esr and esr 
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Since the operator 1Al coincides with that of the environment total spin £„, it turns out that the matrix elements 
(S u r u \0 1Al \ S' U T' U ) arc diagonal 



(S„r u |0 1Al | S' U T' U ) = ^S U (S U + 1)(2S U + l)6 Su ,s>Jr u ,r> u . 
As a consequence, one may check that this property holds for the Hamiltonian of the total cluster 



H 



s' u T' u - l --s'M'vy u ) = 



T(_-[\S' u +S+\/2 / S u 1/2 S 
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A direct calculation shows that the ground state belongs to the Hilbert space sector with S = 3/2 and F = A\. 
Hence, only the environment state with ST = lAi is needed to find the ground state energy (see Table I). 

Let us now consider the next step, an expansion of the current environment block due to the next coordination 
sphere of radius \/2. After an addition of four spins S a , S b , S c , S d , the cluster becomes a square of size 3x3 with the 
bipartite environment of the central site (Fig.l). The basis associated with the added part is 
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(15) 



Repeating the basic steps in the approach we obtain the symmetry adapted basis \SiMj; Tim). The matrix of 
corresponding similarity transformation has the form 



|00;Ail) |00;Sil) \1M;B 2 1) \1M;E1) \1M;E2) \2M-A x l) 

|00: 00) ^ \ 

|11; 00) -\ & 
|01;1M) 0^0^0 
|10;1M) 0^0-^0 

|11;1M) 1 

|11;2M) 1 

The environment Hamiltonian includes only interactions between the first and second coordination spheres 



H n — J 



S ai (3d + Sa^j + (§a + Shj + S 7l (^S b + S c ^j + S m (S c + S d 



(16) 



We now introduce the cluster irreducible tensors W q J^ and U q J^ transforming according to representations T/j, of the 
point symmetry group D4 (for details see Appendix D) 



U q i 1 — — J= (S aq + S b q + S cq + S dq ) , U ql 2 — — — (S aq S bq + S cq S dq ) , 
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and then rewrite Eg. ([IB]) as 
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The reduced matrix elements of the irreducible operators that appear in Eq.([5]) can be obtained exactly from the 
result (fH 
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To compute matrix elements of the Hamiltonian H u we construct the basis 
\ijSjT! ijiSuT u ; S u M u T u pL u ) = 



rat, ran M,MI 



Si Sij S u 
mi mn M u 



Fi Tn r„ 



\iiSimiTi^i) {iiiSiimnTufMn) 



(20) 

formed from the eigenstates \iiSimiTifii) and \inSnmnT nnn) of the "new" and "old" added parts, correspond- 
ingly. Then we obtain using Eci. (f29|) the expression similar to Eq.© with E^SjVi — EinSuFu — 0. Applying exact 
diagonalization to the Hamiltonian H u one can then find the cigcnfunctions 

\i u S u T u fJ, u ) = a i"srr?:i rrSrrTrr l^I^I^I ill^IlT 'u\ S u M u T u fl u ) 

and the energy spectrum £, s r of the environment. By using the recursion relation (for details see (fM|) in Appendix 
B) 

(i u S u r u \\0 1Al \\i' u S' a T' a ) = d ru!T>u ^ E E ^k^iiS^r^^fs^^s^v^ 

ii ,Si ,r i'jjjS'ujF'jj iii,Sn,Tti 



(-l) 1+Sl+S '" +S "[Su,S' u ] 1/2 \ S q ; 1 ^A^iiSiiTiiWO^W^iiS'iiT^) 



(21) 
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TABLE II: Environment states of symmetry QA\ 

iuOA ± iuOAx i u OA x p 
l J± QAi OA) a lElE a 2Ai 2Aj -tW„OAi 

1 0.071 0.449 0.890 -3.651 J 

2 0.569 0.715 -0.406 —0.726 J 

3 -0.819 0.535 -0.205 0.377 J 



TABLE III: Environment states of symmetry 1^4i 



i u lA\ iulAx i u i -A 1 i u 1 A 1 p 

a 0A^lA 3 Q 1A 4 0A 4 a iElE a 2A, 2A, &i u lAi 



1 0.153 0.153 -0.478 -0.851 -3.128 J 

2 -0.470 -0.470 0.566 -0.487 -1.202 J 

3 -0.505 -0.505 -0.672 0.196 1.330 J 

4 -0.707 0.707 



one finds the reduced matrix elements in the environment basis \i u S u A x ) that come into the matrix of the total cluster 
©■ 

The formulas (|3l4l5p allow us to obtain any of possible 54 square cluster states. Our calculation shows that the 
ground state belongs to the Hilbert space sector with S =1/2 and T = A x . Hence, only the environment states with 
ST = 0^2,1^! are needed for the evaluation of the ground state energy. Below we summarize the results obtained 
for this particular case. 

Using Eq.([5]) and the explicit expressions for the nonzero sums of Clebsch-Gordan coefficients of the point group 
Di (see Eq.([27]l in the Appendix A) 

F(A X A X A 1 ;A X ,A X ,A 1 ) = F(EEA X ;A X A X E) = F(EEA X ; EEA X ) = 1, F(A X A X A X ; EEE) = 1/2, 

we obtain 

-73 J 

L / _J 

V3 J J V 3 

-\JlJ -3 J 

in the basis of the states |(L4i0Ai; 0U4i) , \IEIE; 0U4i) , \2Ai2A x ; 0(L4i). The diagonalization of H u { ° Al) yields 
three states of the symmetry (see Table II) As for the 7j„-operator with S u = 1, we have the following matrix 

representation, with the same considerations as for the Hu ^-operator, 

J 
J 

J J 2" J 2 *^ 

-^J -|J _ 

in the basis |0BilBi; IMA X ) , \IB 2 QB 2 ] IMA X ) , |LELE; IMA X ) , \2A X 2A X ; IMA X ). The states of 1A X symmetry are 
listed in Table III. By using the recursion relation (I21|) with the starting value (fl9|) . one finds the reduced matrix 
elements in the environment basis \i u S u A x ). Plugging them into Eq.© we get the target states JiiMA!) [see Eq.©] 
of the cluster and their energies E^^ (i — 1..7). The number of states involved in determining the cluster ground 
state equals 7 (see Table IV). 

Now we list the results for observables. The energy per bond found with the help of Eq.(j4j is e g = —0.3442 J. 
This result may be compared to those results of QMC 7 e g = — 0.3347J, and DMRG e g = — 0.32679J for lattice of 
size 20 x 20 and for number of DMRG states 150&. (Extropolation of the DMRG results in the infinite-lattice limit 
yields e g = -0.3321J). The best available DCM^i, CCM2& and RSRG-EIi2 results are -0.33486 J, -0.33308 J, and 
—0.33409 J, respectively. Using (J6j) we get the ground-state expectation value of the z component of the central spin 

(S'o)o = 0.173 and the staggered magnetization M = i/3(5q)q = 0.299. For comparison, the extrapolated QMC 



n, 



(0A 1 ) 
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TABLE IV: Data on the ground state of the cluster 3x3. 





> J 2 OA, 


Pi OA, 


/3l 9 !A, 


021A, 






E 3 


-0.712 


0.044 


0.010 


-0.695 


0.0048 


-0.090 


0.011 


-4.749 



result for the lattice magnetization M = 0.3070. We also provide an estimate of the spin-spin correlation functions 



(S z S z r=1 ) = -0.115, {S^S z r= ^) = 0.073. 

These estimates should be compared with the known results -0.1116 and 0.0637, correspondingly,—. 

We have made a preliminary calculations by using the small cluster 3x3 and one can see that an accuracy of the 
results is still insufficient. However, we have established the following important features: 

(i) The ground state of the system belongs to identity representation A\ . 

(ii) The lowest-lying environment states of the same point symmetry give a contribution to the ground state of the 
system with 

the largest weight |/3f Ail 2 + I/^iiaJ 2 ~ 0.989. One can see the coefficients /3 2 by nothing that the diagonal matrix 
elements of the reduced density matrix in DMRG language^. 

(iii) A comparison of the ground state energy per bond as calculated by E g /12 = — 0.396J and its infinite-lattice 
approximation ((4]) within our approach shows that we produce a better result. 

At further step, the procedure is repeated and the environment block grows by adding the coordination sphere of 
radius 2. When the new spins S a2 , Sp 2 , Sy 2 , S ri2 of the sphere are added, the cluster transforms into the rhombus of 
size vl3 x \/l3. The cluster has the non-bipartite environment, hence, it is instructive to study this case to examine 
the effect of non-biparticity. 

The Hamiltonian of the new environment decomposes as 



H u — H u (0) + J (^S ai S a2 + SfaSfo + S 71 S 72 + S m S V2 \ . (22) 



H u (0) contains all interactions within the "old" environment, and the second term describes all couplings between 
this part and the added sites. 

The irreducible tensors built from the added spins are the same as those of the first coordination sphere p7|) 



W ql 1 — ^y- {S a2 q + Sj3 2q + S l2q + S V2q ) , W ql 1 — (S a2q S ' p 2q + Sy 2 q S V2q ) , 

W ql — {S a2 q — S l2q ) , W q ® — (S, l2q ~ Sp 2q ) . (23) 

One can then cast the Hamiltonian (1221) in a more amenable form 



Hu = H u (0) + \j [U^ x W 1A ^ + \j [U 1B ^ x H^C 1 + i=J [U 1E x W'Xt 



where U lj are given by 



U ql 1 — ^— (S aiq + Sp iq + S~[ iq + S mq ) , U ql 1 — ^_ (S aiq Sp iq + S liq S mq ) . 



U q l — (Saiq Syiq) i U q2 — (S niq Sj3 iq ) . (24) 

The matrices formed from the reduced matrix elements of W 11 tensor coincide with (fT9|) . To find those of U 11 
tensor we use Eq. (|3~4"]) . The expressions mentioned p^|) are used to initialize the calculations. 

From direct calculations one can show that the quantum numbers S = 5/2 and T = A\ are attached to the ground 
state of the rhombus. This state is formed from 41 envronment states with the symmetry ST U = 2A X and 22 states 
of symmetry ST U = 3A\. Numerical diagonalization gives the cluster ground state energy E g (^Ai) = — 5.779J that 
yields the ground-state energy per bond e g = —0.30925 J in the thermodynamic limit. If we compare this result with 
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TABLE V: Convergence of the ground state properties vs number of environment states kept. 



M(OAi) 


M{lAi) 


E /J 


e/J 


in 


(S§S*(1)) 


(SgS z (V2)) 


1 


1 


-7.9010 


-0.2410 


0.354897 


-0.080333 


0.065647 


5 


5 


-8.1018 


-0.3136 


304148 


-0 1 04533 


071420 


10 


10 


-8.1282 


-0.3238 


0.304928 


-0.107933 


0.073201 


20 


20 


-8.1378 


-0.3279 


0.305707 


-0.109300 


0.073872 


50 


50 


-8.1425 


-0.3301 


0.305101 


-0.110033 


0.074247 


100 


100 


-8.1429 


-0.3303 


0.305187 


-0.110100 


0.074289 


194 


194 


-8.1430 


-0.3304 


0.305187 


-0.110133 


0.074300 



that of QMC, we see that the agreemnet becomes worse. Nevertheless, the conclusions made for the square cluster 
3x3 hold: (i) both the ground state of the environment and that of the total cluster have the lattice point symmetry 
A\. (ii) The largest weight (is of the order 0.993) into the sum of diagonal elements in the density matrix comes from 
three lowest-lying 2Ai states and one state of symmetry 3-Aj, whereas the total number of states is 63. 

Monitoring energies per bond Eisr for the total cluster spectrum Eigp, we found that the minimal value 
£ m in ~ —0.3229 J is reached for the lowest state of symmetry however, E{^Ai) > E g = E{^Ai). A simi- 

lar situation, when a minimal energy per bond belongs to a higher lying state, has been early observed in DMRG 
study of antiferromagnetic chains 2 . Despite the number of sites in the cluster yTS x %/l3 is greater than that of in 
the cluster 3 x 3, we see that the result for e min deteriorates compared to the QMC value —0.3347 J. Close inspection 
allows us to suggest that this is because we are working on the cluster with a non-bipartite environment. 

To proceed with increasing cluster size and satisfy the biparticity requirement we should take the square cluster 
5 x 5 in the next step. For the 24-site environment of the cluster, an exact-diagonalization calculation of the total 
spectrum is not possible at present and so, to move on to the next-larger system, we have to elaborate a procedure for 
determining the states giving the best approximation to true environment states. To solve the problem and implement 
the condition of bipartite environment we take a system in the form of "decorated cross" obtained from the former 
cluster VT3 x \/T3~ by adding four spins S a3 , Sp 3 , S l3 , S V3 (Fig. 1). The form makes equal a number of sites in both 
sublattices, though it incorporates 8 sites that are being attached to the cluster by single lattice bonds. At the same 
time, exact diagonalization of the cluster\/T7 x VT7 is allowed, hence we compare the exact diagonalization results 
with those obtained from a symmetry based truncation procedure and ananlyze a truncation error on a number of 
states kept. Since the cluster increasing is similar to that used in the previous step, we present only the results of 
calculations. The ground state of the extended cluster environment has the symmetry 0A±. The total number of states 
with the same symmetry is 194. Together with 439 lAi-states of the environment they form a ground state of the 
total cluster labeled by the symmetry numbers \A\. Results for the ground state energy per bond e = —0.3304, the 

staggered magnetization m = 0.305 and the spin-spin correlation functions (SqS^ =1 ) = —0.1101, ^o^^.-V2 / = 0-0615 
agree well with the mentioned ED and QMC results and are much better than those obtained for the square cluster 
3 x 3. A deviation from the ED result is found for (SqS* =2 ) = 0.0169. This discrepancy arises from finite size effects 
and an imperfect topology of the cluster. 

We now describe the low-energy spectrum of the environment. As the dynamics of Neel order parameter is the one 
of a free rotator, the low-energy levels scale as E(S) ~ S(S + 1)/N, where the inertia of that rotator is proportional 
to the number of site a 28 i 29 . The environment lowest-energy levels (tower of states) belonging to different irreducible 
representations of the lattice point group are shown in Fig. 2 for different S sectors. The SU{2) breaking due to long- 
range Neel order appears as a set of Ai-states, lying off from other levels, with an energy scaling as E(S) ~ S(S + 1). 

In the remainder of this section we describe a version of the truncation procedure. The main idea will be illustrated 
on an example of the ground states properties. An inspection of results for the current and previous clusters reveals 
that one have to take the lowest-lying environment eigenstates both in the 0Ai and lAi sectors. As for the number 
of kept states it seems to be most simple to take M states equally from the both subspaces, albeit the choice may 
not be optimal. To prove that this concept works we recalculate the observables found above on various number of 
envronment states kept (see Table V). As can be seen from Fig. 3 the convergence of the results is exponentially fast 
in M. Merely keeping 100 basis states may be as efficient as keeping of all 633 environment states intact. We regard 
the resulting better than 0.01% agreement for e and m as support for the efficiency of our truncation procedure. 
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IV. CONCLUSIONS. 



In this paper we present a quantum cluster solver for spin-S" Hcisenberg model on a two-dimensional lattice. The 
formalism is based on the real-space renormalization procedure and uses the lattice point group-theoretical analysis 
and nonabelian SU (2) spin symmetry technique. Let us summarize advantages of the approach: 

(i) The cluster spin states are decomposed into parts belonging to different irreducible representations of the lattice 
point group and to different values of the total spin. Due to the embedded group-theoretical analysis, our approach 
can handle each of the cluster target states independently that offers a distinct advantage for parallel computation. 

(ii) An extension of MPM destined for quantum spin chains to higher dimensions has inspired construction of 
variational methods for the ground states of 2D spin Hamiltonians (vertex state models^, tensor product variational 
approach 35 , tensor product ansata 3 ^). Since, the trial states are represented by two-dimensional product of local 
weights, these approaches are face with severe limitations concerning their applicability because of relation between 
a spin value and lattice topology. The shortcoming lacks in our formalism. 

(iii) Large sparce-matrix diagonalization algorithms (Lanscoz technique, for example) used in DMRG and ED 
methods converge to maximum and minimum eigenvalues of a model Hamiltonian, i.e. to eigenvalues at the edges of 
the spectrum. Our approach gives access to eigenstates of an entire spectrum. 

(iv) Combined with decimation procedure of the environment states like those used in DMRG technique the group- 
theoretical analysis allows us to overcome exponential growth of computational efforts with increase of system size. 
Our approach using the total spin S and the irrep index T/i as good quantum numbers yields a rather reliable 
truncation procedure of the Hilbert space of the model Hamiltonian. 

(v) Calculation of observables for the central spin involving a density matrix of the environment reduces edge effects 
which are inevitable on finite-size clusters. 

The major drawback of the formalism is that it does not allow an easy implementation: a complexity in construction 
of basic sets via repeated evaluation of 6j and 6r symbols, the calculation involves two matrix diagonalizations etc. 
The performance gains from implementing the SU(2) and lattice point symmetries are not impressive in comparison 
with gains from exploiting just the simple U(l) symmetry leading to total magnetization as good quantum number. 
Their using in studies with larger clusters without truncation cannot help to alleviate the problem of exponential 
growth of computational efforts. 

In the method that we suggest, short range correlations on the scale of the cluster are taken into account, while 
correlations on a scale larger than the cluster size are neglected. To overcome this shortcoming we need to restore 
translational symmetry of the lattice. The results of these investigations will be reported elsewhere. In this connection, 
we note that the translational invariance holds for the DCM and CCM methods. 
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V. APPENDIX A 



Let \iiSiTi iiiSuTii; SMTfj,) is a state with total spin 5, third component M, and transforming according to 
irreducible representation Ffx. This state appears in the tensor product decomposition (iiSjTi) x (inSnTn), where 
(iST) denotes a state with total spin S, irreducible representation T and i labels other possible quantum numbers. 

We need to compute the matrix element 



iiSiTi ijjSijTjj; SMTfi 



E E 

qv {m.fi} 



7 7 M 








r 


* 




r' 






' Si 


Sii 


S 




' S'j 


D n 


S" 


V v 1 








M _ 












mi 


run 


M 






m'jj 


M 



{iiSjmiTiHz \ Ugj\ ijSimiVjfi'j) \iiiS u muTuHu 



IF 



i 7 
qv 



(25) 
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The Wigner-Eckart theorem for a double irreducible tensor reads 



(i I S I m I T m \W^\i\S' I m\T' I ^) = (-l) s ^l Sl 1 S \ ) 

V —mi q thj I 

where the 3j symbol is related to the Clebsch-Gordan coefficient by 



7 r' 7 r, 



(i I s I r I \\w 1 i\\i' I s' I r' I ), (26) 



Si S2 S3 



1 



=(-1) 



Si+S 2 -m 3 



Si S2 5*3 
772 x 777,2 — m-3 



I 777i rn 2 7773 I V2S , 3 + 1 
A full contraction of five Clebsch-Gordan coefficients of the point group may be written via the 6r symbol 



F(r J r /J r ; r , J r / Jj7 )= £ 



7 7 Ai 




Ti 


Tn 


r 


* 




r' 


r " 




"7 r. 


Ti 


* - 


V v 1 




. fl1 


Hii 


m _ 






M77 








P-i _ 





7 r^z r„ 

f // /7 fin 



Ti 


T/7 


'} 


"a 


r' 


r 


r' 
1 77 


r'7 




1 


/<' 


/<■ _ 



(27) 



however, it is more convenient to find directly this sum. 

Substituting (12"6")) into ([2!>]) and performing the sum with the aid of Eq. (|2"?| and the formula (see^, for example) 

x ^r \^ -a X)\-X P rj 



we get finally 



r a s 
—p a' a 



s p t 
-a —ip —t 



(-!) a I q P a 
(2a +1) I s r t 



iiSiTi inSiiTu; SMF/j, 



[if 1 -* x W 1 ^ 1 i^S'jT'j I'uS'jjT'jj-S'M'T'ii' 



(28) 



dss'SMM'Srr'Spii' (— 1) 



Sjj 1 



F(r J r JI r ; r' J r'„7) 



x (ijSjTj \\U^\\ i'jS'jT'j) {mSnTn \\W^\\ i' JI S' II T' tI ) 
The reduced matrix elements appearing in (|29[) result from the previous iteration. 



(29) 



VI. APPENDIX B 



The systematic increasing cluster size requires an iterative procedure to compute the reduced matrix elements of 
the double irreducible tensors C/ 17 or W 1 " 1 (acting on the states with indices I and //, respectively) in the basis 



Si 


Sii 


S 




Ti 


T/7 


r 


m 7 


777// 


777 






M77 





li/Siro/IV/jr) \iiiSiimnTnfiu) 



(30) 



with aid of the Wigner-Eckart theorem. On the other hand one can use the basis of states (|3"0|) to obtain 

(iSmFn\W£}\i' 'S'm'V '//) 
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a i I s I r i; i II s II r II u -t I s I r I ;i' II s' II r' II 



J' ST' 



E 



mimum T 



Si Sh S 
mi rrin m 



Si Sjj S 
mi m'jj mf 



(-1) 



Sii—mii 



SlI 1 Sjj 

~mn q m'jj 



ii T// r 

fl! [in P 



ri t'jj r 
[II fJt'u n 



7 r' J7 r w 



The sum over m/, to// and to^j is performed with the aid of the formula 



(iuSuTuWW^Wi'uS'nT'jj) 



EH-) 

X</>P 



,« -V+?-x+r-p [pa r \ ( r c p 

Kip a -x J \ X P ~P J \ P 7 



(31) 



a b c 
—a —(3 —7 



a b c 
r p q 



(32) 



The sum of three Clebsch-Gordan coefficients of the lattice point group in turn can be transformed as follows 

Tl T'jj F 



E, E 



Ti r u r 

pi pn p 



pi Pn p 



7 r' 7/ r 7/ 

v h'jj mi 



17 r /7 r 

A* 



E 



7 r r 

V fl' p 



E # E 



Tl T'jj r' 

M// A' 



17 r JZ r 

At/ A*// A 4 



7 r w 
i/ ^' 7/ mi 

r, r JZ r' 

At/ A*// A*' 



E 



7 r r 

v [i! p 



7 r 7/ r// 

" A*// A*// 



- * 


7 


r' 


r 




i> 




Al _ 


- * 


7 


r 


f 






p' 


p 



After permutation of the first and second columns in the third Clebsch-Gordan coefficient the sum over projections 
At/, fin, fijj, and is easily performed that gives immediately 61" symbol 38 



E 

re 



" 7 r 


r 


'i 


v p! 


m _ 





7 r' r// 



r\ r r 



;/ ^^^(r/r^r'), 



(33) 



where we use the symmetry property of the Clebsch-Gordan coefficients 

= e(rir 2 r) 



ri r 2 r 

p\ pi p 



r 2 r 1 r 

P2 Pi P 



and the sign e (ril^T) — ±1 depends on the point group. 

The reduced matrix element can be computed using (|31I33|) that yields the results 



x (-i) 1 + s *+ s i*+ s [5,5'] 



S 1 S' 



i' ST' 

;iiiS II T II a i I S I T I ;i'S'T> 



1/2 



(34) 



and 



(iSTWU^Wi'S'T') 



/ j fx iiSiV 1 \iuS u Vn ,x %' l S' I V' l \%uSuVii 



x(-l) 1+Sl+s " +s '[S,S'} 1/2 \ S , 1 S ' \(i I S I T J \\U 1 ni' I S' I T f I )\ 7 ^ ^ 



(35) 



where [S] = (25+1). 
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VII. APPENDIX C 



Here, we give a detailed derivation of representation for the D4 point group in the basis of four coupled connector 
spins 



Si 52 S12 
mi m 2 m 12 



ISMS^SaS^S^SM) = 
Consider first the ir/2 rotation C4 about the z axis 
C , 4|5i52(5i 2 )5 3 5 4 (534)5M) = E E E 



S3 S4 S34 

m 3 7774 TO34 



mi,m2 m3,m4 77112,77734 



5*1 S*2 S*12 

mi ui2 m\2 



S12 S34 S 
mi2 m 34 M 



5*3 5 4 S34 

TO3 TO4 m34 



|5imi) |5 2 to 2 ) \S3m3) \S4m4) 



S\2 S34 S 

m\2 m-34 M 



\S2m2) \S3m3) \S4m4) \Si 



= E E E 

mi, 7712 7773,7774 77712,77734 



Si S*2 S*12 

mi m 2 mi2 



S*3 S*4 S34 

W3 m4 m34 



5i2 S34 S 
mi2 m 34 M 



x ^ ^ |S 2 m 2 ) |S 3 m 3 ) |5 4 m 4 ) |5imi) 

^m 2 777 2 ^777 3 777 3 ^777 4 777 4 ^777 i 7 

777 1 ,7772 7773,7774 

To proceed we use the orthogonality property for the Clebsch-Gordan coefficients 



^777 2 7712^77737773^77747,. : 



^77717771 ^ ] ^ ] 

Sl4, 77714 S23, 77723 



52 ^3 523 

m 2 m 3 m 23 



52 53 523 

m 2 m 3 m 2 3 



54 5i 5i4 
m 4 mi m 14 



54 5i 5i4 
7714 fh\ m\4 



Then 



C 4 15x52(5x2)5354(534)5^/) = ^ Yj E 



mi,777 2 7773,7774 77712,77734 



5x 52 5x2 
mi m 2 mi2 



53 54 534 
m 3 m 4 m34 



5i2 534 5 
TO12 m, 34 M 



x E E l 5 2W2) |5 3 m 3 ) |5 4 m 4 ) |5imi) V 



7771,7772 7773,7774 



S14, 77714 S231 m 23 



52 53 523 
7772 »773 m23 



52 53 523 
77l 2 m 3 77723 



54 5i 5i4 

TO4 777 1 777 X4 



54 Si 5x4 
TO4 777x mi4 



- E E E 

7771,7772 7773,7774 77712,77734 



S\ S2 5x2 

777 x 777 2 m,X2 



53 54 534 
7773 7774 77734 



5x2 534 5 
777x2 m, 34 M 



Y E I* 527 " 2 ) \ S 3 m 3) \S4m4) |5x777x) 



777i, 7772 7773, 777 4 



x E E E 

S14, 77714 523,77723 77714,77723 



52 53 523 
777 2 "7 3 77723 



52 53 523 

777 2 m 3 m 2 3 



54 Si 5x4 

TO 4 777 x 777X4 



54 Si 5x4 
7774 777x m 14 



^777 2 3 777 2 3 ^777 1 4 777 1 4 



E E E 

7771,7772 7773,7774 77712,77734 



5x 52 5x2 
mi m 2 mx2 



53 54 534 
7773 7774 777.34 



5x2 534 5 
777x2 m 34 M 



Y E l^ 2 ^ 2 > l^ms) |5 4 777 4 ) |5xmx) 

mi ,m2 m3,m4 



x E E E 

S14, 77714 S23, 77723 77714,77723 



52 53 523 
777 2 f«3 W23 



52 53 523 
777 2 m 3 77723 



54 Si 5x4 

7774 777x mx4 



54 Si 5x4 

7774 777 x 771x4 



E 

S,M L 



523 5x4 5 
77723 mx4 M 



523 5x4 5 
77723 TO14 M 
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A full contraction of six Clebsch-Gordan coefficients yields the 9j symbol 

E E E 



mi,m2 rri3,m4 mi2,m34 _ 



Si 5*2 <Sl2 

mi m 2 mi2 



5*3 5*4 ^34 
TO 3 TO4 7T134 



5*12 S34 S 1 
77J12 m 34 M 



5*4 Si S14 
TO 4 TO! TO 14 



5*2 S3 S23 

to 2 m 3 m 23 



S23 5*14 S 

Tl23 Wl4 M 



^ ^^5l+54 — 5l4 ^ + — S34 ^ ^ 



S23+S14 — S 



E E E 

mi,m2 ui3,m4 mi2, 77134 



Si 54 5i4 
mi m 4 mi4 



5*2 5*3 5*23 

m 2 m 3 m 23 



Si S2 S12 
mi m 2 m i2 

5*14 S23 S 
TO14 m 2 3 M 



S3 S4 S34 

TO 3 TO4 TO34 



5*12 S34 5* 

TO12 m 34 M 



{Si S2 S12 

Si S3 S34 

S14 5*23 S 

where we use the symmetry relation for the Clebsch-Gordan coefficients 

_ ^_-^"jSi + S , 2 — S12 

By introducing the state with the recoupled four spins 

|S 2 S3(S23)S 4 Si(Si4)5M) = E E E 



Si S2 S12 
mi m 2 m i2 



52 Si S12 
m 2 mi m i2 



mi,m2 ni3 mi4,m23 



5*2 5*3 S*23 
TO 2 TO3 TO23 



S4 Si S14 

m 4 mi mi4 



523 5i4 S 

m 2 3 wi4 M 



x |S 2 m 2 ) |S 3 to 3 ) |S 4 to 4 ) |Simi) , 



we get finally 



C^ISi^S^Sa^C^SM) = E (-lf 1+Si - Sl4 (-lf 3+S4 - S34 (-lf 23+Sl4 -s 



514-523 



x [(25i2 + l)(25i4 + l)(25 2 3 + l)(25 3 4 + l)] 1/2 { S 4 S 3 S 34 

Sl4 S23 S 



Si S2 S12 



} \S 2 S 3 (S 23 )S 4 Si(Si 4 )SM) , 



that is transformation from a coupling scheme to another. Then an action of the operator C4 is defined by a linear 
transformation of the basis k' = { S[ S' 2 ( S' 12 ) S' 3 S' 4 ( S 34 )} 

C i \„)=Y J D^ K {c i ) |«'), 

k' 

where the matrix D^ K (Cij is determined by the expression 

D^l (4) - Ss^Ms'Ms'Ms^s'Js^su [(2S12 + 1) (25i4 + 1) (2S 23 + 1) (2S 3 4 + I)] 1 



/2 



^ ^ 2^5l+54 — S14 ^ -j^S^ + S^ — S34 ^ ^S23 + S'l4 — S 

We can handle analogously another symmetry operations in the D4 group. 



Si S2 S12 

5*4 S'3 S34 

5*14 S23 S 
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VIII. APPENDIX D 



The standard method for constructing an irreducible basis is to use the projection operator 



[ J sec 

and the shift operator 

[ J g£G 

where [G] is the order of the group G, [T] is the dimension of the irreducible representation T and Dj^J (g) are the 
irreducible matrix elements, fj, or v is an index enumerating the basis. Supposing that ip is one of the reducible basis 
vectors of G, an irreducible basis might be obtained by applying 

^=«-^)<- (38) 

If is the basis for the irrerp T then 

Plu€ = ^l- (39) 

Let us construct irreducible tensors from the operators {S a , Sb, S c , Sd} forming the nearest environment of the 
central site. The transformations of one of the given spins under the elements of the group are 

E S a = S a , C4 S a = Sd , C4 S a — S c , C4 S a = Sf, , 

G% S a — Sb j C2 S a = Sd , C v S a = S a , C v S a = S c , 
that together with (|36p gives immediately the irreducible tensors of the one-dimensional representations 

U\t = N Al (S a + S b + S c + S d ) , U\t = 0, 



Uqi 1 = 0, = N B2 (S a -S b + S c - S d ) . 

To find the irredicible basis for the two-dimensional representation E we construct the projection operator Pf and 
then apply it to the spin S a that yields 

Uql = Ne (S a + Sb — S c — Sd) ■ 
Using the shift operator ([57]) and acting according to the rule (|59")l we obtain the second irrep basis vector 

^q2 = Ne [S a — Sb — S c + Sd) ■ 

We choose the coefficents Nr (T = Ax, E) so that Hamiltonian written through the irreducible tensors coincides with 
the initial spin operator form. 

It is not always possible to construct all irredicible basises from one chain. The theory says that one have to 
choose another starting function. In a computer realization, therefore, we build all chains generated by all vectors of 
a reducible basis 

.,/,-, ^ /).„.:,/:,•,. (40) 

3 

and form the matrix X T from the chains 

^i = NrY. D l^9)D ji {g). (41) 
9 
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The rank of this matrix 

1 ' g 

determines a number of linear independent columns, where the character x (<?) = (ff)- After these columns are 

established, orthogonalized and normalized with the help of Schmidt-Gram procedure we get the first Cr columns of 
the transformation matrix 2$ rv- By running over all irreducible representations and repeating the basic steps in the 
approach we obtain the square matrix of corresponding similarity transformation onto the symmetry adapted basis 

i 

As an example we calculate ^5^g 34 f° r ^ ne nearest-neighbor environment of the central site. The characters of 
3-dimensional representation can be read off from the 3x3 matrices given in Sec. Ill We thus obtain 

E C4, , Cf Cf Cf , C| cr^ , 
X 3 -1 -1 -1 1 

whence we conclude = D^ 1E \ A direct calculation of A r matrices from Eq. pTj) yields 

/ i i 



2 2 

I I 
\0 





1 


1 


\/2 




4 

1 


4 

1 


4 


I 


4 
n/2 

4 


4 

V2 
4 


4 

1 

2 



and 




By noticing that ranks of the matrices equal to unity, we find finally via the Schmidt-Gram procedure the transfor- 



mation matrix T, 



(i) 



\IM-Bx) \\M\E1) \1M\E2) 



01; 1M) 


1 


1 

2 


1 

2 


10; IM) 


1 


1 
2 


1 
2 


11; 1M) 





1 


1 

V2 



IX. APPENDIX E 

For the reader convenience we give the character table of the group D<± 





E 


c^c- 1 


c l 


°2 i u 2 




Ax 


1 


1 


1 


1 


1 


A 2 


1 


1 


1 


-1 


-1 


Bx 


1 


-1 


1 


1 


-1 


B 2 


1 


-1 


1 


-1 


1 


E 


2 





-2 
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and the matrices of double irreducible representation taken in the basis xy (see22, for example) 
^ (£) (Cf ) = ( J _°i ) • D(E) ^ = ( ^ 1 ) ' = ( 1 ) ' = ( -1 ^ ) ■ 
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Fig.l Clusters used in the calculations. 

Fig. 2 The lowest-energy spectrum of the environment for the cluster \/l7 x y/l7 on the square lattice. The SU(2) 
symmetry breaks and a long-range Neel order appears as a set of ^i-states with an energy scaling as E(S) <~ S(S + 1) 
(dashed line). The symbols represent the irreducible representations of the different eigenstates. 

Fig. 3 The cluster ground state energy E, the energy per bond e, and the staggered magnetization m convergence 
for the y/Tf x y/Vf cluster vs number of environment states kept. 



